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Abstract,  Models  of  a  coupled  bilinear  structure  are  studied  here.  Included  are  appli¬ 
cations  to  immunology  and  generalized  defense  processes.  In  some  cases,  coupled  bi¬ 
linear  systems  arise  naturally  and  in  others  they  represent  valid  approximations  to  more 
highly  nonlinear  systems.  A  basis  is  developed  for  multi-level  command  and  control. 
Finally,  hierarchical  pursuit-evasion  strategies  are  introduced,  and  plans  are  presented 
for  more  detailed  studies. 

Keywords.  Bilinear  control;  cybernetics;  game  theory;  hierarchical  systems;  and 
model  1  ing. 

INTRODUCTION 


The  purpose  of  this  paper  is  to  analyze  the  hier¬ 
archical  structure  of  general ized  defense  and 
pursuit-evasion  games  which  are  amenable  to 
bilinear  or  linear  system  analysis.  Also,  coupled 
bilinear  systems,  C8LS,  sometimes  arise  in  a 
natural  manner  for  complex  physical  processes.  In 
other  cases,  they  may  be  developed  to  approximate 
input-output  dynamics  of  a  more  highly  non!  inear 
nature.  Cellular  fission,  nuclear  fission,  convec¬ 
tive  and  conductive  heat  transfer  and  certain  chem¬ 
ical  reactions  are  good  examples  of  the  former, 
while  rigid-body  dynamics,  such  as  appear  in 
aircraft,  ship  and  robotic  control,  are  good 
examples  of  the  latter.  Even  for  many  of  these 
cases,  however,  the  CBLS  representation  may  be 
exact  by  a  redef inition,  and  an  increase  in  the 
number,  of  state  variables.  BLS  and  CBLS  models 
for  some  of  these  are  studied  by  Mohler  (1973). 

Time  and  space  limitations  only  permit  a  prel im- 
inary  analysis  of  this  complicated  problem  along 
with  a  systematic  overview.  However,  a  base  is 
established  for  future  research  with  particular 
application  to  generalized  defense  processes  (human 
and  national). 


STRUCTURAL  PROPERTIES 

For  this  paper,  BLS  are  given  by  the  state 

equation, 

Ilf  =  Ax  +  £  BkV  +  Cu,  (1) 

k»l 

where  x  CRn,  u  £  Rm  with  components  uj ,  ..., 

Un,,  and  A,  C,  ,  k  *  1,  ...,  m,  are  appro¬ 

priate  real  constant  matrices.  Sometimes,  it  is 
convenient  to  define  the  output  of  the  BLS  (1)  by 


y  »  Dx,  (2) 

where  y  C  RP,  and  0  is  an  appropriate  real  con- 
stant  matrix. 

It  is  readily  seen  that  addending  the  state  vector 
x,  with  xn+1  =  1,  allows  (1)  to  be  written  in 

the  homogeneous  form  (HBLS);  i.e.,  (1)  with 
zero  C. 

Then  it  is  shown  by  Brockett  (1973)  that  systems 
with  state  polynomials  (products  of  states)  and 
quotients  of  states  may  be  decomposed  into  new 
HBLS.  Consequently,  it  is  noted  by  Tenney  (1984) 
that  systems,  which  involve  transcendental  func¬ 

tions  of  the  state,  may  be  decomposed  into  "poly¬ 
nomial  systems"  by  merely  defining  new  state 
variables  as  the  functions  themselves.  This 

generalization  of  BLS  makes  them  and  certain  of 
their  associated  methodologies  amenable  to  a 
larger  class  of  nonl inear  systems  which  includes 
rigid-body  dynamics  such  as  for  maneuvering  air¬ 
craft,  ships,  land  vehicles  and  robots.  Also, 

predator-prey  processes,  combat  games,  irrmunology, 
el ectrical-power  generation,  and  heating-ventilat¬ 
ing  and  air-conditioning  (HVAC)  processes  have 
similar  representations.  AH  of  these  systems  are 
quite  complex  with  hierarchical  structures  of 
appropriate  orders. 

Within  each  BLS  itself  there  evolves  a  certain 
hierarchical  structure  which  is  associated  with  a 
canonical"  decomposition  such  as  convenient  for 
its  Vol terra  representation  (assuming  its  exist¬ 
ence,  of  course).  In  this  manner,  BLS  (1)  and  (2) 
may  be  generated  by 

y.Dx.ogx,*^*,,  (3) 

1*1  i=l 

where 
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•  m 

*2  ’  Ax2  +  E  Bkukxl  +  Cu> 

k*l 

(4) 

x\  =  Axi  +  E  VkVi  +Cu> 

k=l 


with  x(0)  =  0  for  convenience. 

It  is  readily  seen  [Bruni,  et  al  (1974);  Mohler 

(1973)]  that  Xj,  1  =  1,  ...,  correspond  to  the 

terms  In  the  Vol terra  series,  and  the  correspond ing 
kernels  are  generated  as  a  “nesting"  [Rao  and 
Mohler  (1975)]  of  1  Inear-system  impulse  responses 
according  to  (4). 

The  Vol terra  series  for  BLS  may  sometimes  be 
approximated  by  a  finite  number  of  terms  (4),  and. 
In  the  case  of  weakly  BLS,  may  be  exactly  repre¬ 
sented  by  a  finite  nunber,  or  decomposed  Into  a 

finite  nunber  of  linear  systems  with  outputs 
multiplied  together  to  form  Inputs  to  successive 

1  Inear  systems  according  to  (4),  Consequently,  a 
physical  system  approximated  (or  given)  by  a  finite 
hierarchical  structure  associated  with  (4)  is 
conveniently  analyzed  by  1  inear-system  theory. 


PHYSICAL  PROCESSES 

The  hierarchical  surveill ance,  decision  and  control 
structure  which  evolves  for  CBLS  is  demonstrated 
here  by  several  relevant  examples.  Figure  1 
Indicates  the  natural  model  which  Is  derived  from 
the  chemical  law  of  mass  action  and  from  cellular 
fission  for  the  immune  system.  See  Mohler,  Brunl 
and  Gandol  f  1  (1980) ,  and  Mohler  (1982). 


immunity  and  a  future  secondary  response  to 
Infection. 

Hierarchical  control  even  occurs  In  the  generation 
of  Ab  classes.  E.g.,  one  antibody  is  more  effec¬ 
tive  early  In  an  immune  reaction  while  another  Is 
effective  later.  Hence,  the  well-known  switch-over 
phenomenon  for  IgM  to  IgG  Is  shown  to  be  somewhat 
optimum  by  Perel son  (1977).  Anatomically,  there  is 
a  hierarchical  preference  in  the  circulation  of 
Immune  function.  B  and  T  cell  precursors  migrate 
to  the  lungs  rapidly,  but  they  also  spend  only  a 
short  duration  there.  On  the  other  hand,  a  large 
number  migrate  to  the  spleen  more  slowly,  and  are 
trapped  there  for  a  much  longer  time  period. 

Many  of  the  immune  control  mechanisms,  which  are 
associated  with  the  nonlinear  feedback  block  in 
Fig.  1,  are  not  well  understood.  But  the  plant 
Itself  is  pretty  well  defined  by  CBLS  as  shown  In 
Fig.  1.  Detailed  mathematical  models  are  presented 
by  Mohler,  Brunl  and  Gandol  fl  (1980),  and  Mohler 
(1982). 

Just  as  the  Immune  response  evolves  in  a  hierarch¬ 
ical  fashion,  so  does  the  evolution  of  the  human 
from  the  fertilized  egg.  A  single  cell  is  pro¬ 
grammed  genetically  with  the  necessary  prescrip¬ 
tion,  and  life  evolves  by  a  sequence  of  carefully 
controlled  chemical  reactions,  cellular  divisions 
and  differentiations.  Again,  CBLS  formulate  the 
basis  of  a  mathematical  synthesis. 

National  defense  takes  on  a  hierarchical  C^ 
structure  similar  to  that  of  body  immune  defense. 
Naval,  air  and  land  forces  are  coordinated  with 
each  also  having  subforces  to  deal  with  alien 
intrusion.  Central  command  headquarters  may  al¬ 
locate  forces  and  various  defense  mechanisms 
according  to  certain  criteria  of  performance. 

At  the  final  level,  competition  between  defense 
forces,  Xj  and  yj,  similar  to  biological 
species,  may  be  modeled  by  CBLS  in  the  general  form 


First,  the  human-defense  system  includes  various 
mechanisms  to  detect  and  to  destroy  an  alien  which 
is  similar  in  structure  to  that  associated  with 
national  defense.  These  include  the  humoral  system 
(generating  antibody,  Ab),  the  cell-mediated  irmune 
system  (CMI,  generating  macrophages,  M^,  cytotoxic 
T  cells,  T r,  helper  T  cells,  T^,  and  suppressor 
T  cells,  TSJ,  and  stages  of  complement. 

The  latter  consists  of  nine  cascaded  BLS  (i.e., 
enzyme-protein  reactions)  which  successively 
amplify  enzyme  concentration  until  it  "drills"  the 
alien  (e.g.,  tumor  cell).  Rreaking  down  the  alien 
permeab 11 ity,  fluid  leaks  In  until  the  alien  cell 
may  burst  apart. 

Ab  receptors  on  the  surface  of  certain  white-blood 
cells  detect  aliens  of  specific  chemical  structure, 
and  generate  appropriate  defense  mechanisms  to  best 
isolate  and  destroy  the  alien  substance  (e.g., 

bacteria,  chemical  contaminants ,  tumors,  etc.). 
Certain  mechanisms  are  ineffective  against  certain 
aliens.  E.g.,  Ab-complexes  (Ab,  and  MJ  may 
even  shield  a  tumor  while  complement  (triggered  by 
a  specific  Ab)  may  be  effective.  The  system 

involves  chemically  specific  and  nonspecific  sur¬ 
veillance,  communication,  command  and  control 
(C3).  M<£,  Tc  and  Ab  seem  to  be  capable  of 

alien  destruction  on  their  own,  but  the  most 
effective  defense  involves  an  orchestrated 
combination  of  these  processes. 

Feedback  control  evolves  according  to  alien  concen¬ 
tration  and  chemical  structure,  and  according  to 

concentrations  of  Ab,  Th  and  Ts.  The  latter 
turn-on  and  turn-off  the  humoral  process  according 
to  its  need  for  Ab  and/or  memory  cells  to  establish 


Tt  -  -V  +  Wx  +  V 

$  ■  -Ayx  ♦  By(uy)y  ♦  y 


(5) 


Here  x,  y  represent  opposing  forces  of  appropriate 


composition  and  dimension; 


are  appro¬ 


priate  vector  force  replenishments;  Ax,  Ay  are 
appropriate  matrix,  opposlng-attr ition  coefficients 
such  as  for  "aimed  fire";  Bx(-),  By(-)  are  ap¬ 
propriate  matrix,  self-attrition  and/or  generation 
coefficients.  ux,  u»  are  appropriate  cross 

controllers.  For  "area  Tire"  of  constant  level  per 
force;  Bx(*),  By(* )  are  diagonal  matrices,  such 
that  the  1th  diagonal  element  is  the  ith  control 
component.  I.e.,  uxj  *  -axj  -  dxy. 
Here,  d^  Is  an  appropriate  constan' 
and  axj  is  a  self-attrition 

de 


I.e 

ux  Is  an  appropriate  constant  row  vector, 
axj  is  a  self-attrition  coefficient. 
Uy|  is  defined  similarly.  This  CBLS  is  a  gen¬ 


eralized  representation  of  the  traditional 
Lanchester  model.  Dolansky  (1964),  and  Wozencraft 
and  Moose  (1984). 


Hierarchical  CBLS  appear  In  numerous  engineering 
applications  such  as  chemical  processes,  nuclear 
power  generation  and  heating  ventilating  and  air- 
conditioning  (HVAC)  systems,  see  Mohler  (1972).  In 
the  latter,  for  example,  bilinear  heat-transfer 
models  are  prominent.  Decisions  are  made  to  use 
solar  storage,  heat  pumps  and  evaporative  or 
refrlgeratlve  air  conditioning  according  to  changes 
in  weather  and  various  energy  demands  and  costs. 
The  latter  Involves  stochastic  processes  as  do  most 
of  those  discussed  above.  For  HVAC  systems,  random 
weather  predictions  are  significant  in  the  command 
and  control  process,  and  actual  random  weather 


variations  can  affect  the  model.  Similarly,  obser¬ 
vations  of  alien  state  Is  statistical  in  nature. 
The  typical  sonar  or  radar  tracking  in  national 
defense  and  chemical  tracking  in  immune  defense  are 
excellent  examples.  The  synthesis  of  CBLS  filters 
for  this  problem  is  studied  by  Halawani,  Mohler  and 
Kolodziej  (1984),  and  is  shown  to  offer  better 
performance  than  the  extended  Kalman  filter. 

While  stochastic  games  are  significant  for  such 
applications,  only  the  deterministic  case  is  con¬ 
sidered  next,  as  a  prelude  to  the  even  more  diffi¬ 
cult  problem.  Here  competition  may  take  place 
between  defense  forces  and  alien  forces. 


PURSUIT-EVASION  GAMES 

Several -pursuers,  one-evader  deterministic  games 
are  the  subject  of  this  section.  The  solution  to 
the  game  requires  the  Introduction  of  a  strategic 
variable,  then  the  game  of  kind  is  studied  to 
determine  possible  cooperations  between  two  pur¬ 
suers.  The  optimal  solution  to  a  linear  guadratic 
game  and  a  simpler  form  of  this  solution  are  given, 
yielding  the  general  solution  to  the  N  Pursuer  - 
one  evader  game.  A  sensitivity  study  leads  to  a 
simple  hierarchical  structure  which  greatly  reduces 
the  amount  of  computation  required. 

let  the  position  (or  more  generally  state),  Xj, 
of  an  evader,  E  with  control  v,  relative  to  pur¬ 
suer,  Pi  with  control  Uj,  be  given  by 


Here,  Xj,  Uj,  v  and  fi  have  appropriate 
dimensions. 

The  capture  occurs  whenever  one  pursuer  Pj  gets  E 
within  its  terminal  manifold  expressed  by  the  quad¬ 
ratic  equation  at  terminal  time,  t fi 

XVMTxlW  -  r1’ 

where  Txj  is  a  given  matrix  and  ri  a  scalar. 
The  team  of  pursuers  minimizes  an  appropriate  per¬ 
formance  index  J,  which  E  maximizes.  This  defines 
a  single  Hamiltonian  for  the  team  of  pursuers. 
Therefore  whenever  a  pursuer  is  added  to  the  game, 
one  equation  is  missing.  The  corresponding  missing 
variable  will  be  called  the  strategic  variable 
Zj,  for  player  Pj.  This  variable  will  be 
defined  such  that  when  Zj  *  0,  Pj  is  of  no  use 
to  the  solution  of  the  game,  and  the  larger  Zj, 
the  more  effective  Pj  is. 


Ill)  z3pj0(p j/pj*) •  Pj0  belongs  to  this  zone 

then  Pj  will  play  a  role  in  the  (Pj,  Pj,  E) 

game. 

z1Uz2Uz3  =  zPjo(Pj  >.  Pj)  =  the  case  in  which  the 

outcome  of  the  (Pj,  E)  game  is  lower  than  for 
the  (Pj,  E)  game. 

1v)  z4Pj0(pj/pi*) ■  If  pj0  belon9s  t0  this  zone 

then  Pj  will  play  a  role  in  the  (Pj,  Pj,  E) 

game. 

v)  ZgPjo^j/Pj*) •  If  pj0  belongs  to  this  zone 

then  Pj  can  play  a  role  in  a  game  involving 
Pj  only  if  E  does  not  play  optimally  or  if  the 

presence  of  the  pursuers  forces  E  to  deviate  from 
the  (Pj,  E)  optimal  strategy. 

vi)  z6Pj0(Pj/P j*) •  If  PjQ  belongs  to  this  zone 

then  Pj  will  not  play  any  role  In  any  game  in¬ 

volving  Pj,  provided  that  Pj  plays  optimally. 
In  this  zone,  Zj  =  0. 

z^Uz^Uzg  *  zPjo(Pi  Pj)  *  the  case  in  which  the 

outcome  of  the  (Pj,  E)  game  is  lower  than  for 

the  (Pj,  E)  game.  The  six  zones  cover  the  whole 
space.  Here,  Pj,  Pj  and  Pj  refer  to  three 

different  classes  of  copursuers  in  the  game. 

As  an  example,  the  N-pursuers,  1-evader  version  of 
the  wall-pursuit  game  presented  by  Isaacs  (1965) 
could  be  renamed  N  cutters  and  fugitive  railroad, 
where  the  evader  of  maximum  speed  1  is  constrained 
to  the  y  axis.  The  minimum-time  problem  is  con¬ 
sidered,  where  the  pursuers  Pj  and  ?2 ,  con¬ 
trolling  their  constant-speed  heading,  have  the 

following  respective  characteristics:  the  capture 
sets  are  circles  of  radii  1  and  0.5;  the  maximum 
speeds  are  2  and  1.5,  and  Pin  *  (-4,  4).  See 
Fig.  2. 

The  zones  are  the  same  for  every  pursuer  identical 
to  Pj  but  are  a  function  of  They  were 

derived  from  the  solution  to  the  1  vs.  1  game  and 
according  to  geometric  considerations.  the  same 
kind  of  definitions  can  be  made  for  3  pursuers; 
e.g., 

‘i'WV'Y’  pj*> 

etc.,  but  the  derivation  of  the  zones  would  require 
that  the  2  vs.  1  game  be  solved  first. 


The  Game  of  Degree 


The  Game  of  Kind 

Once  the  setting  of  the  game  is  known,  those  pur¬ 
suers  that  will  be  relevant  must  be  determined. 
One  way  to  proceed  is  to  define  for  each  pursuer 
the  area  where  a  copursuer  must  be  located  in  order 
to  allow  (or  not)  cooperation.  For  each  pair  Pj, 
Pj  the  space  is  divided  into  the  following  six 
zones: 

1)  zipj0(p"j/pj*^  •  pj0  belongs  to  this  zone 

then  Pj  will  not  play  any  role  in  any  game  in¬ 
volving  Pj,  provided  that  Pj  plays  optimally. 
In  this  zone,  Zj  =  0. 

ii)  z2PjQ(P^/Pj*).  If  PjQ  belongs  to  this  zone 

then  Pj  can  play  a  role  in  a  game  involving 
Pj  only  if  E  does  not  play  optimally  or  if  the 
presence  of  other  pursuers  forces  E  to  deviate  from 
the  (Pj,  E)  optimal  strategy. 


The  game  studied  is  a  linear,  quadratic,  general¬ 
ized  2-pursuer,  1 -evader  game  with 


dxi 

dt  *  A1x1  +  8 iu1  +  C1v’ 


(7) 


where 


Xl(tf )TTx1*i{tf)  *  r?,  and 


\  j  <2ZxiQixi +  u{Riut  -  vTsv)dt 


1 


(8) 


is  to  be  minimized  by  the  pursuers  and  maximized  by 
the  evader. 


1-pursuer  vs.  1-evader  game. 


Let  Xpj,  \ej 


be  the  costate  vectors  required  by  the  game,  such 
that  the  transversal ity  condition  is  given  by 

xpi(V  "  z1Tx1xl(V*  *  zl'Tx1x1*tf' 

where  z^,  z/  can  be  any  nonzero  positive 
scalars.  The  Hamiltonians  are 


+  CjS  C2Te2<#>2(t,  t0) 

.  (13) 

*2(t,  tQ)  *  (Az  +  N2Tp2  +  H2Te2)<i2(t,  tQ) 

+  C2S  ■c|Tel(®1(t ,  tQ) 


Hp  =  (A1X1  +  B1ui  +  C1v)  xpi 

-  "2  ( x  .jX  ^  ”  x  3v ) , 

He  *  (A1X1  +  Biu1  +  C1v)Txei 

+  j  (x^Q1x1  +  -  vTSv) . 

The  open-loop  Nash  solution  to  the  game  Is 

V  ■  Vxio  > 

v*  =  S  tQ)x1o  , 


(9) 


(10) 


V**  V  “  <A1  +  N1Tp1  +M1Te1^1(t-  V' 


where  <^(t,t)  =  I 


N1  =  =■  C^C},  and 


•lfT 


tpi  -*A;Tpi  -  Vi  -  vitpi  -  TPiMiTei +  °i- 


(id 


Tei  *  -AiTei  '  Te1A1  *  Te1N1Tp1  *  Te1M1Te1  ’  <V 
where 

WV  *  z1Tx1*  WV  *  zi'Txr 

To  this  solution,  another  two  conditions,  restrict¬ 
ing  the  solution  must  be  added.  I.e.,  z-j  >_  0 
when  substituted  into  the  Hamiltonian  provides  the 
playability  condition,  and 


Xp  i  ^  tf ) x 1 ( tf )  1  0 
is  the  capture  condition. 


The  solution  is  a  classical  one.  A  complete  deri¬ 
vation  of  the  two-player,  linear,  quadratic  game  Is 
given  by  Ichikawa  (1976)  and  Hamalalnen  (1978).  A 
procedure  to  decouple  the  Riccatl  equations  (11)  is 
given  by  Simaan  and  Cruz  (1973)  taking  advantage  of 
a  prel iminary  solution  cornnon  to  several  Initial 
conditions,  reducing  the  problem  to  the  computation 
of  successive  linear  equations. 

2-pursuer  vs.  1-evader  game.  The  four  costate 
vectors  required  are:  Xpj,  Xej,  i  =  1,  2, 

where 


xpi^V  *  z1Tx1x1^V  *  xei*V  =  zl’Tx1xl'V 

The  open-loop  Nash  solution  to  the  game  is 

U1 =  "WptV1*  Vxio«  1  *  *• 2-  (12> 

V*  =  s  (CjTg^jft,  tQ)x10  +  C2Ie2  2(t,  to)x20), 
^(t,  tQ)  =  (Ax  +  N1Tpl  +  MiTel^l^t*  to^ 


where 


^(t.t)  »  *2(t,t)  =1,  and 


TplXl’[-TplAl-Arpl-'TplNlTpl-TplMlTel 
+  Ql]xl  ■  TplClS  lc2Te2x2’ 


Te1xl  *  -  A'T*i  '  T-iNiT"i  *  T^MiT. 


el  1  1  el  'eri'pl  ‘el 'Tel 

-  -  TgjCjS  1C2Te2x2, 


Vx2  “  ^”Tp2A2  *  A2Tp2  '  Tp2N2Tp2  *  Tp2M2Te2 


(14) 


+  Q2]x2  -  Tp2C2S 

Te2x2  *  C-Te2A2  -  A2Te2  -  Te2N2Tp2  -  Tg2M2Te2 
-  Q2]x2  -  Te2C2S 


where 

W  ■  Z1TX1’  W  3  Z*1TX1-  1  *  2- 

The  solution  is  very  complex  compared  to  the  1  vs. 
1  game.  The  differential  systems,  (13)  and  (14) 
are  tightly  coupled  requiring  parallel  computation 
whereas  the  solution  to  the  1  vs.  1  game  could 
proceed  in  several  steps. 


Note  that  only  one  of  the  two  pairs  of  strategic 
variables  (zj,  zj),  ( Z£  >  zz)  can  be  fixed, 
the  other  one  must  be  selected  properly  to  cor¬ 
respond  to  the  initial  condition  of  the  game. 


Keeping  the  same  controls,  but  setting  z?  *  0 
z?  *  0  and  exploring  all  the  possible  xgftf) 
will-  give  the  help  zone,  H  =  Z1UZ2UZ3UZ4, 
the  area  outside  of  which  does  n°t  play  a  role 
in  the  optimal  (Pj,  Pg,  £)  game. 

Simplifying  the  2  vs.  1  solution,  define  such 
that 


Z1Telxl 


Te1x1 


1  *  1,  2. 


(Note  that  Zj  *  I,  the  identity.)  Now,  the  solu¬ 
tion,  while  (12)  and  (13)  remain  identical,  is  as 
follows. 


Tpl  ’  'TplAl  '  AlTpl  ‘  TplNlTpl  '  TplMlTel 
-  TplClS‘lc2Z2Tal  +  V 


Tel  "  -TelAl  *  AlTel  *  TelNlTpl  *  TelMlTel 


(15) 


-  TgjCjS  1C2Z2Tel  -  Qj, 


where 


W  *  zlTxl’  Tel(tP  *  Z1T; 


xl’ 


Tp2  '  'Tp2A2  *  A2Tp2  '  Tp2N2Tp2  "  Tp2M2Te2 


*  Tp2C2S  lciZ2lTe2  +  Q2 


Te2  *  ”Te2A2  '  A2Te2  *  Te2N2Tp2  '  Te2M2Te2 


-  Te2C2S  1c}z21Te2  -  Q2 


where 


(16) 


l.n 


Let  6,  -  -Tf  ^  Zj  ZilTi 

j=1 

1.1-1 

E  ZjZ1T1+W1 

J 

1+1, n  , 

and  Wi  -  -T|M  ^  Zj  Z‘JiT1 

j 

Except  for  the  addition  of  0-s  and  the  sum  In  v*, 
the  solution  Is  Identical  to  tne  1  vs.  1  game. 


Tp2^V  3  z2Tx2 '  Te2^f^  =  z2Tx2;  and 

Z2  -  Z2(a{  +  QjT*j)  -  (A*  +  o2t;|)z2. 

Z'1  =  Z^CaJ  +  02T~|)  -  (a{  +  (17) 

WTxlW  *lfTx2W>  Z21(tfKZ2(tf)]-1. 

The  similarities  with  the  1  vs.  1  game  are  now 
obvious;  the  equations  are  decoupled  and  can  be 
solved  in  turn.  However  the  problem  of  guessing 
the  strategic  variable  remains.  Note  the  simplifi¬ 
cation  that  Qj  *  Q2  3  0  would  introduce. 


This  suggests  a  sensitivity  study  on  both  0 ^  and 
W-j.  In  some  instances  the  0-j  are  negligible 
when  the  pursuers  represent  equal  threats  to  the 
evader.  In  this  case  the  difference  between  the  1 
vs.  1  game  and  the  n  vs.  1  game  is  seen  to  come 
hardly  from  the  controls  adopted  by  the  pursuers 
but  from  the  control  adopted  by  the  evader.  bj 
is  an  important  term  if  P-j  plays  a  minor  role  in 
the  capture.  On  the  other  hand  is  always  neg¬ 
ligible  when  the  playability  condition  is  not 
violated,  provided  that  the  pursuers  be  classified 
in  order  of  importance,  being  the  most  impor¬ 
tant  (the  closest  one  to  E  for  the  minimum  time 
problems).  Consequently,  the  structure  is  shown  in 
Fig.  4  where  each  pursuer  solves  (18)  -  (20), 
finding  the  couple  x«j(tf),  z\  corresponding 
to  x<f0,  and  passes  this  information  forward. 


Consider  Fig.  3  as  an  example  of  a  2 -dimensional , 
linear,  quadratic  game  with  identical  pursuers. 
Here,  Q  =  0,  R  *  S  —I, 


r-i  ii 

'2' 

V 

it 

ii 

CO 

C  = 

r- 

i 

r\> 

i _ 

i — 
o 

1 _ 

-0. 

Also,  uj,  U2»  v  are  scalars,  Tx  *  I,  r  =  1. 

The  capture  zone  for  the  1  vs.  1  game,  defined  by 
J  >_  -lQlO  is  very  limited,  as  is  the  usable 
part  of  the  terminal  manifold.  However,  for  the 
2  vs.  1  game,  there  is  no  limitation  on  the  usable 
part,  and  the  capture  zone  is  much  larger.  Ending 
on  identical  points,  variations  on  the  strategic 
variable  Z2  leads  to  different  trajectories.  The 
proper  values  for  x^ftf),  xgCtf)  and  Z2 
must  be  selected  in  order  to  fit  with  the  initial 
conditions  xjq,  X2o* 

Solution  to  an  N- identical-pursuer,  1-evader  linear 

quadratic  game.  The  solution,  for  0=0,  takes  the 

following  form  with  Xf(t),  Xj(t0)  *  xjg. 
given  by  (7): 


A  parallel  structure  can  be  derived  from  this  "rip¬ 
ple  structure,"  enhancing  the  independence  of  the 
individual  pursuers  with  respect  to  the  team,  but 
at  the  expense  of  an  increased  number  of  equations 
to  solve.  A  "minor"  pursuer  can  be  added  but  the 
gain  produced  by  this  pursuer  must  be  weighted 
against  the  amount  of  delay  or  computation  that 
this  very  pursuer  will  have  to  cope  with.  This 
structure,  implemented  for  the  previous  example, 
did  not  show  any  variation  in  state,  control  or 
performance  superior  to  0.1%  with  the  rigorous 
sol ution. 


CONCLUSIONS  AND  DIRECTIONS 

Convenient  hierarchical  structures  for  pursuit- 
evasion  games  (or  command  and  control)  are  gener¬ 
ated  along  with  an  overview  of  the  role  of  BLS  in 
generalized  defense  and  other  relevant  applica¬ 
tions.  Future  plans  call  for  an  integration  of 
this  analysis  with  stochastic  nonlinear  observa¬ 
tion,  and  stochastic  BLS  models  such  as  are 
relevant  to  estimation  and  control  in  the  applica¬ 
tions  presented  above. 


U1  =  bTt1x1*  v*-C  EVf.  U8) 


*r.  .  -T.A  -  ATT,  -  VH  +  N)T(  -TjM  V  zj  z:1^, 

FI  (19) 
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u  refers  to  control  reactivity  with  subscripts 
c  for  complement,  B  for  B  cells,  T  for  T 
cells,  M  for  macrophages. 

BISC  refers  to  coupled  BLS. 

C3,  Cg  refers  to  appropriate  complement 
stages. 

Other  terms  defined  in  text. 

Fig.  1  BLS  Synthesis  of  Significant  Immune  Proc¬ 
esses 


ralroad 


P20  C  Z1  : 

:  P2  alone  will  capture  E. 

P20  c  h  '■ 

:  ?2  alone  will  capture  an  optimal  E 

P20  C  Z3  : 

:  P2  will  be  helped  by  P^  in  captur- 

ing  E. 

P20  C  Z4  : 

:  P^  will  be  helped  by  P2  in  captur¬ 

ing  E. 

P20  C  Z5  : 

:  Pj  alone  will  capture  an  optimal  E 

P20CZ6  : 

:  P^  alone  will  capture  E. 

Fig.  2  Influence  of  the  Position  of  ?2  on  the 
Resulting  Game.  P^q  =  (4»  -4)« 


Fig.  3  2-Pursuer,  1-Evader  Relative  Trajectories 
Leading  to  Capture. 


Fig.  4  n-Pursuer  Suboptimal  Hierarchical  Struc¬ 
ture. 


